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' Abstract 
<^ . 

, We analyze the horizon structure of famiUes of space times obtained 

by evolving initial data sets containing apparent horizons with several 
connected components. We show that under certain smallness conditions 
' the outermost apparent horizons will also have several connected com- 

, ponents. We further show that, again under a smallness condition, the 

maximal globally hyperbolic development of the many black hole initial 
qh| data constructed in [9], or of hyperboloidal data of [19], will have an event 

horizon, the intersection of which with the initial data hypersurface is not 
^Jjl connected. This justifies the "many black hole" character of those space- 

times. 



^ ' 1 Introduction 



There is an ongoing effort to construct "many black hole" solutions of the vac- 
uum Einstein equations numerically (see e.g. [2, 25, 26] and references therein). 
In practice this means that one numerically evolves initial data which contain 
trapped surfaces for as long as the computer allows. The question then arises 
whether the resulting space-time does indeed contain more than one black hole, 
or for that matter, any. Several issues arise here: 

a) The notion of a black hole is usually tied to the existence of a conformal 
completion of the space-time (but see [8] for alternative proposals). It 
is far from clear that the vacuum solutions, which are in principle asso- 
ciated with their numerical counterparts discussed in [2,25,26], possess 
sufficiently controlled conformal completions, if any. 
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b) Even assuming the issues in point a) do not occur, consider an initial data 
set g, K) which contains several trapped, or marginally trapped, sur- 
faces. If yet another trapped or marginally trapped surface Sq encloses all 
the previous ones, then the geometry enclosed by Sq is hidden from ex- 
ternal observers by the null hypersurface J~^{Sq). Numerical calculations 
of tentative radiation patterns inside J'^(S'o) have absolutely no relevance 
to the data detected at . Thus if one is willing to associate black hole 
regions to apparent horizons, then only the outermost apparent horizons 
are relevant. In this context a condition for a multi-black-hole space- 
time would be that the outermost apparent horizon has more than one 
component. 

c) In any case, the event horizon itself might have nothing to do with the 
apparent horizons, even the outermost ones. Under appropriate hypothe- 
ses, the existence of an apparent horizon implies the existence of a black 
hole region, but this can be much larger than the region enclosed by the 
outermost apparent horizon. In particular one might imagine a situation 
in which the outermost apparent horizon has several components, all of 
which are enclosed by a connected event horizon, so that the space-time 
contains only a single black hole region. 

The object of this paper is to point out that the issues raised above can be 
analysed in a reasonably satisfactory way for the "many Schwarzschild" initial 
data constructed in [9], or for the data obtained by the gluing constructions of 
[19, 21], or for families of initial data sharing certain qualitative properties, as 
made precise below, similar to those of [9, 19]. The first main result here is that 
for a rather general class of "small-data" families of black-hole space-times, 
the outermost apparent horizon will have several connected components. 
We prove this both on the usual asymptotically flat initial data hypersurfaces 
and on hyperboloidal ones. For the initial data of [9] the relevant smallness 
condition holds when the mass parameters mj, z = 1, ... ,1, of the individual 
Schwarzschild black holes are small enough as compared to the distance pa- 
rameters Tj. For the initial data of [19] the smallness condition holds when the 
gluing necks are sufficiently small. 

One of the features of the initial data sets of [9] is that these metrics are 
exactly Schwarzschild outside a compact set, and this guarantees that for any 
one of these, the associated maximal globally hyperbolic development {^,g) 
necessarily has a which is complete to the past.^ As already indicated, the 
existence of J^~^ is the usual starting point for a discussion of black hole regions. 
The second main result here is the proof that, for certain configurations and 
again for mass parameters small enough, the intersection 

^+ := j-(j^+)n^ (1.1) 

^We note that both here, and in several situations of interest, one can use the results in [22] 
to infer the existence and past-completeness of J^"*". However, the estimates there are based 

on spherical outgoing null hypersurfaces, which can be used to prove the existence of, at best, 
a connected black hole region (if any). Further, the differentiability properties of J^^ which 
can be directly inferred from that work (compare [7]) are not sufficient to be able to invoke 
the stability results of [15], as needed below. 
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of the future event horizon J'(,y+) with the initial data hypersurface has 
at least / components. (Indeed, we show that si has exactly I components and 
believe that this should also be true for (f^; a proof of such a claim about S"^ 
would require complete control of the global structure of the resulting space- 
time, which is well beyond the range of techniques available nowadays.) 



2 "Many black hole" initial data 

There are several constructions of families of initial data containing apparent 
horizons, see [1,3,4,14] and references therein. In this section we briefly de- 
scribe three such families of "many black hole initial data" . Before doing this, 
it is useful to recall how apparent horizons are detected using initial data (com- 
pare [4]): let, thus, {<9^,g,K) be an initial data set, and let 5 C ^ be a 
compact embedded two-dimensional two-sided submanifold in S^. If is the 
field of outer normals to S and H is the outer mean extrinsic curvature^ of 
S within ^ then, in a convenient normalisation, the divergence 9^ of future 
directed null geodesies normal to S is given by 

0+ = H + Kij{g'^ - n'n^) . (2.1) 

In the time-symmetric case 9+ reduces thus to H, and S is trapped if and only 
if i7 < 0, marginally trapped if and only ii H = 0. In the hyperboloidal case 
with Kij = Qij we obtain 9^ = H + 2. 



rrii 



2.1 Brill-Lindquist initial data 

Probably the simplest examples are the time-symmetric initial data of Brill and 
Lindquist. Here the space-metric at time t = takes the form 

g = ^l;^{dx^ + dy^ + dz^) , (2.2) 

with 

I 

1=1 

The positions of the poles Xi G and the values of the mass parameters G R 
are arbitrary. If all the rrii are positive and sufficiently small, then there exists 
a small minimal surface with the topology of a sphere which encloses Xi. This 
follows from general arguments in geometric measure theory, as implemented 
and described in more detail in Section 3 below. In addition, from [22], the 
associated maximal globally hyperbolic development possesses a ^+ which is 
complete to the past, but the differentiability properties of the conformally 
completed metric may not be sufficient to justify some key steps of the global 
analysis below concerning event horizons. 



^We use the definition that gives H = 2/r for round spheres of radius r in three dimensional 
Eudidean space. 
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2.2 The "many Schwarzschild" initial data of [9] 

There is a well-known special case of (2.2), which is the space-part of the 
Schwarzschild metric centred at xq with mass m : 

V 2|x-xo|/ 

where S is the Euclidean metric. Abusing terminology in a standard way, we 
call (2.3) simply the Schwarzschild metric. Allowing the mass parameter to 
be nonpositive leads to naked singularities or flat regions, cf. [9], but we shall 
always require positive masses here. The sphere |x — xo| = m/2 is minimal, and 
the region \x — xo\ < m/2 corresponds to the second asymptotic region, beyond 
the Einstein-Rosen bridge. 

Now fix the radii < 4i?i < i?2 < oo. Denoting by B{d,R) the open 
coordinate ball centred at a with radius R, choose points 



Xi e ro(4i?i,i?2) := 



B{0,R2)\B{0,ARi) , Ri>0 
B{0,R2), Ri = 



and radii r^, i = 1, . . . , 2N, so that the closed balls B{xi, 4ri) are all contained 
in ro(4i2i,i?2) and are pairwise disjoint. Set 

n := To{Ri,R2) \ (uiB{xi,ri)^ . (2.4) 

We assume that the Xi and arc chosen so that is invariant with respect to the 
reflection x — > —x. Now consider a collection of nonnegative mass parameters, 
arranged into a vector as 

M = {m, mo, nil, ■ ■ ■ ,m2N), 

where < 2mj < Vi, i > 1, and in addition with 2mQ < i?i if i?i > but 
mo = if -Ri = 0. We assume that the mass parameters associated to the 
points Xi and — are the same. The remaining entry m is explained below. 

Given this data, it follows from the work of [13] (as pointed out in [9], 
compare [10]) that there exists a (5 > such that if 

2N 

y\mi\<6, (2.5) 



i=0 



then there exists a number 



2N 



m = ^mi + 0{S^ 



and a C°° metric which is a solution of the time-symmetric vacuum con- 
straint equation 

R{9m) = , 

such that: 
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1. On the punctured balls B{xi,2ri) \ {xj}, i > 1, ^^g is the Schwarzschild 
metric, centred at x^, with mass m^; 

2. On \ 5(0, 2i?2), gg^ agrees with the Schwarzschild metric centred at 0, 
with mass m; 

3. If Ri > 0, then agrees on B{0,2Ri) \ {0} with the Schwarzschild 
metric centred at 0, with mass mo. 

In fact, this construction also gives that is symmetric under the parity map 
X — > —X. 

2.3 Black holes and gluing methods 

A recent alternate technique for gluing initial data sets is given in [19], see 
also [21] for the time-symmetric case and [20] for more general results in the 
asymptotically Euclidean case. In this approach, general initial data sets on 
compact manifolds or with asymptotically Euclidean or hyperboloidal ends 
are glued together to produce solutions of the constraint equations on the 
connected sum manifolds. Only very mild restrictions on the original initial 
data are needed. The neck regions produced by this construction are again of 
Schwarzschild type. The overall strategy of the construction is similar to that 
used by Corvino (and in many previous gluing constructions). Namely, one 
takes a family of approximate solutions to the constraint equations and then 
attempts to perturb the members of this family to exact solutions. There is a 
parameter rj which measures the size of the neck, or gluing region; the main 
difRculty is caused by the tension between the competing demands that the 
approximate solutions become more nearly exact as r/ — > while the underlying 
geometry and analysis become more singular. In this approach, the conformal 
method of solving the constraints is used, and the solution involves a conformal 
factor which is exponentially close to 1 (as a function of rj) away from the neck 
region, but which is nonetheless not completely localised. 

Consider first an asymptotically flat time-symmetric initial data set, to 
which several other time-symmetric initial data sets have been glued by this 
method. If one assumes that the resulting necks are mean outer convex, as de- 
scribed in detail in Section 3 below^, then the existence of a non-trivial minimal 
surface, hence of an apparent horizon, follows by standard results, cf. Section 3. 
This implies the existence of a (possibly disconnected) black hole region in the 
maximal globally hyperbolic development of the data. As for the Brill-Lindquist 
construction, the asymptotically flat initial data produced in this way may not 
have sufficient differentiability at the resulting to obtain good information 
about S'^. 

Consider, next, hyperboloidal initial data with 

Kij = gij . (2.6) 
^This will hold if the gluing regions are made small enough. 
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It follows from (2.1) that in this setting trapped or marginally trapped surfaces 
are characterised by the condition 

e+ = H + 2<0 . (2.7) 

Fixing a polar coordinate r on the standard three-dimensional hyperboloid, the 
constant curvature —1 metric takes the form 

g = -T^dr^ + r^dn"^ 
+ 1 

(where dO.'^ is the constant curvature +1 metric on S"^), then it is straightforward 
to calculate that the mean curvature with respect to the outer normal of the 
'constant r' geodesic spheres is given by the formula 

= 2Vl + r-2 . 

Now suppose we glue together two hyperboloidal initial data sets. From the 
point of view of far away observers sitting on the other side of the ensuing neck, 
the inner pointing normal for a geodesic sphere on one half of this configura- 
tion is actually pointing towards them, thus outer-pointing as far as they are 
concerned; hence the quantity 

-H + 2 = -2/(r^ + + 1) 

measures "trapedness" with respect to the other asymptotic region. This is neg- 
ative for any r > on the hyperboloid, and will remain strictly negative when 
r is large enough even after the gluing has been performed. This means that 
large spheres on one side of the neck are trapped from the point of view of the 
^+ on the other side and hence, by standard Lorentzian geometry arguments, 
can not be seen from that Scri. This again implies the existence of a black 
hole region. Notice that this simplest case is rotationally invariant, and there 
is a unique minimal sphere encircling the neck (sec Lemma 3.4 below). Hence 
by continuity, there is at least one marginally trapped, i.e. with 6+ = 0, rota- 
tionally symmetric geodesic sphere. A similar argument establishes existence 
of black hole regions when several initial data sets satisfying (2.6), at least two 
of which are hyperboloidal, are glued together. 

3 Outermost apparent horizons with many compo- 
nents 

In this section we generalise the examples above and consider a family of asymp- 
totically Euclidean metrics {grj}, < V — VOi which satisfy the two properties 
below. Our goal is to show that when ijq is sufficiently small, the outermost 
apparent horizon of each of these metrics has a large number of components. 

We assume that for r) G (0,770], {^,gr)) is a Riemannian manifold with 
boundary of dimension 3 with a single asymptotically Euclidean end E, and 
such that is a union of / copies of 5^. (In the context of the initial data 
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of Section 2.2, this amounts to removing from the manifold that part which 
hes on the other side of the connecting necks.) We suppose furthermore that 
around each boundary component there is an annular 'neck region' Aj^, i = 
equipped with a diffeomorphism : x [—1,1] ^ Ai, such that 
^i{S^ X {-1}) = (dy) n Ai. Thus 

y = E{r]) U ^1 U . . . U ^7 

is a union of manifolds with boundary, intersecting only along the submanifolds 
{dE{r})) r\ Ai = ^i{S'^ X {+1}) so that the Ai are mutually disjoint. We call 
$1(5*^ X {-!}) and $j(S'^ x {+!}) the outer and inner boundaries of Ai. The 
end E[rj) is diffeomorphic to an exterior region in R^, and we fix a family of 
diffeomorphisms 

: \ UiB{xi,piiri)) E{ri), 

and assume that the radii of these balls Pi{r)) tend to as 77 \ 0. (It is only 
a matter of convention that we think of the annular regions as fixed, whereas 
E{r]) is identified with an r/-dependent region. However, the metric gjj varies 
nontrivially on each of these regions.) 

Our hypotheses on the metrics g"^ are as follows: 

a) [Metric convergence on the distinguished end:] If K is any compact 
subset of \ Ui{xi}i=i^... J, then for some a G (0, 1) 

}i^l\\Ki9n) - ^Wc^'^^iK) = 0; 

here 6 is the Euclidean metric on M^. 

b) [Mean outer convex necks and small minimising cycles:] For rj in 

a sufficiently small interval (0,170), both the inner and outer boundaries 
(5^ X {±1}) of Ai are mean outer convex with respect to g^j] further- 
more, there exists a smoothly embedded sphere Si which represents the 
fundamental class (jj G -^2(^1, Z) and with area \Si\ ^ as rj ^ 0. 

Each of the three constructions outlined in Section 2.2 produce families 
of metrics satisfying these hypotheses. For example, for the construction in 
Section 2.2, if Mq := {mQ,mi, . . . ,m2N) is a (2A^ + l)-tuple of nonnegative 
numbers and M{r]) = (771(77), r/Mg) is the associated mass-parameter vector from 
that construction, then g^j := ^^g^^^ satisfies both these hypotheses. Similarly, 
the initial data of Section 2.3 satisfy the hypotheses here if 77 is a sufficiently 
small parameter controlling the outer radii of the I necks across which the 
gluing is performed. 

We begin with a geometric result which holds under slightly more general 

hypotheses: 

Lemma 3.1 Let g be a Riemannian metric on A = S'^ x [—1, 1] such that the two 
boundaries S'^ x {±1} are mean outer convex. Fix a generator oa for H2{A,Z). 
Then any surface S which is absolutely area minimising in this homology class is 
smoothly embedded, lies in the interior of A, and consists of a single component 
of multiplicity one. 
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Proof: The existence of a homological area-minimiser S in the class of integral 
currents in a manifold with mean outer convex boundaries, and the regularity of 
its support, is a standard result in geometric measure theory, cf. [30, Theorems 
37.2 and 37.7]. (These arguments work equally well for domains with mean 
outer convex boundaries, cf. [29], and by the maximum principle, the support 
of the resulting minimiser is disjoint from dA.) In particular, the support of E 
is a finite union of smooth, oriented, connected surfaces Si, ... , Ej, where each 
Sj appears with some non- vanishing integer multiplicity kj . Thus on the level 
of homology 

ki[^i] + . . . + kj[^j] = aA, 

whereas 

|E| = |fci||Si| + ... + |fcj||Ej|. (3.1) 

We claim that the support of S has only one component, and this occurs 
with multiplicity 1. To prove this, note first that any component divides 
S'^ X [—1,1] into precisely two components. This may be seen by 'capping off' 
the boundary S*^ x {— 1} of ^ by adding a 3-ball; the interior of the resulting 
manifold AuB^ is diff'eomorphic to M^. By the Jordan separation theorem, any 
smooth, oriented, connected surface Ej embedded in A, hence in M^, divides 
this space into an 'inside' and an 'outside'. For example, a point p lies in the 
inner component if (all) generic paths 7 connecting p to the outer boundary 
5^ X {1} intersect Ej an odd number of times. In any case, this decomposition 
shows that in homology, [Ej] = ±aA or else [Ej] = for each j. If any Ej 
is null-homologous, then we can obviously discard it, since it adds a positive 
amount to the area of E without contributing to the homology class; possibly 
changing orientations, we can therefore assume that each [Ej] = aA- 

Finally, amongst the Ej select one, E', with smallest area. Then from (3.1), 
|E'| < |E|, and equality holds only if E' is the only component, and occurs 
with multiplicity 1. Thus E' is the connected homological area-minimiser, as 
required. □ 

Now let us return to the more general situation. Using this lemma, we 
represent the generator aj = [^j{S'^ x {+1})] of H2{Aj,Z) by a homologically 
area minimising surface E^; according to hypothesis b), aj is also represented 
by the sphere Sj. Both Ej and Sj are smoothly embedded, connected surfaces 
of multiplicity one. (Since grj has nonncgative scalar curvature, it is known 
[5, 29] that Ej - or indeed any stable minimal surface - must be either a sphere 
S'^, or possibly a torus if grj is flat in a neighborhood of Ej.) By assumption, 
\Sj\ 0, and hence |Ej| as well. 

It is proved in [18, Lemma 4.1] that with the hypotheses above, for every 
< 77 < ryo there exists a unique outermost minimal surface S^, which is a 
union of embedded stable minimal spheres of class C^'+i'" jf g,^^ jg of class C^'". 
Furthermore, if we denote by the exterior of Sjj in y (i.e. the unbounded 
component of \ S^), then 5^ is absolutely area minimising in its homology 
class in and moreover, is simply connected. 
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Theorem 3.2 There exists rji e (0, ?7o] such that if r) e (0,771], then 5^ C 
uf^-^^Ai and the intersection of 5^ with each annular region A-i is nonempty. 
Hence has at least I connected components. If we assume that there do 
not exist any stable minimal homologically trivial surfaces in any of the regions 
{Ai,gjj) when 77 is small enough, then S^HA^ contains exactly one component, 
and hence S^j has precisely I components. 

Proof: Let £'(0, R) denote a large sphere in which contains all of the points 
Xi, and let Q denote the part of ,y interior to this sphere. Coherently orienting 
the fundamental classes aj{H2{Aj,Z)), we have that [^(O, i?)] = a := ai + 
... + (77, where we regard aj G H2{Aj,Z) ^ H2{S^,'L), as induced by the 
inclusions Aj ^ S". Prom [18, Lemma 4.1], we know that S^' is diffeomorphic 
to the complement of a finite number of spheres in M"', and hence S*^ must 
be homologous to 5(0, R) as well, i.e. [5^] = a. For each ry we choose area- 
minimising representatives Sj (r/) of Uj in Aj , as in the preceding Lemma. By 
hypothesis a), 5'(0, R) is mean outer convex for 5^ if 77 is small enough, since it 
is strictly convex for the limiting Euclidean metric S. Thus we have 

IS1I + ... + IS7I < |5J < \S{0,R)\. 

The first inequality holds because USj is absolutely area minimising in its ho- 
mology class in ,5^, while the second inequality follows from the fact that S*^ 
is absolutely minimising in its homology class in J^'. We claim that for rj suf- 
ficiently small, >S',j lies in the union yli U . . . U ^7. Granting this claim for the 
moment, let us prove that 5^ has at least / components. Choose for each j a 
smooth embedded curve 7^ which connects the inner boundary $j(5'^ x { — 1}) of 
Aj to 5(0, i?), does not intersect any of the other annular regions Ai, i ^ j, and 
which represents the Poincare dual of aj in Hi[Q,d^). Then the homological 
intersection number of 7^ with [5^] equals 

([7j-],^) = ([7j],c^i + --- + ^/) = l- 

On the other hand, if ■~fj is in general position, then this intersection number is 
also computed by counting the signed geometric intersections of this curve and 
this surface. Therefore this geometric intersection is nontrivial, which shows 
that Srf n 7^ for each j, and hence Srf has at least / components. 

To prove the claim, suppose there exists a sequence r]i ^ such that S{i) := 
contains a point G O \ UAj with ^ ^ G \ {xi, . . . , X7}. The interior 
curvature estimate for embedded stable minimal surfaces proved by Schoen [28] 
states that there is a uniform upper bound for norm squared of the second 
fundamental of S{() with respect to g^^ near q^^. More precisely, for any p G S{t) 
with p{p) = minj{||? — pi(r])\} > 5 > Q ioi £ sufficiently large, there exists a 
constant C > 0, independent of £, such that \IIs((.){p)\'^ ^ C. By standard 
calculus, this imphes that the portion of S{t) in a ball of radius pip/l) around 
p may be written as a graph with uniformly bounded gradient over a disk of 
radius p{p)/^ in T^S{1). In particular, the area of S{tj is uniformly bounded 
below by a positive constant. 

Applying these bounds to a finite covering of Q\\JiB{xi, p) for any p > 0, and 
then taking a diagonal subsequence for some sequence pj — 0, we may extract a 
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subsequence S{i') which converges to a nontrivial smoothly embedded minimal 
surface S{oo) in M'^ \ {xi, . . . Since all of the S{i') are unions of spheres, 

and the number of components is uniformly bounded, the limiting surface must 
have finite genus. In addition, S{oo) is compact and has bounded area. We may 
now apply a well-known removable singularities theorem for minimal surfaces, 
see [6, Prop. 1] for a proof, which shows that S{oo) is a nontrivial compact 
embedded minimal surface in R^. Since no such surfaces exist, we have reached 
a contradiction. We have now proved the first assertion, and hence that 5^ has 
at least one connected component in each Ai. 

For the remaining assertion, write Si{T]) = Srj H Ai, and suppose that this 
surface has more than one component for some i, i.e. Si{r]) = Uj^iSij{r]), where 
J > 1 and the Sij{ri) are smooth embedded surfaces. By the same argument as 
in Lemma 3.1, each Sij{r]) separates A^ into two components. If A^ contains no 
null-homologous stable minimal surfaces, then each component of Ai \ Sij{ri) 
must contain exactly one of the two boundaries ^i{S^ x {±1}). However, the 
components Sij{r]) are disjoint, and so if there are at least two, then any one 
must be contained in either the interior or exterior region of another; since 
their union is an outermost surface this is impossible. We conclude that Si{ri) 
is connected. This completes the proof. □ 

In the case of data of Section 2.2 the hypotheses of the second part of 
Theorem 3.2 are verified: 

Corollary 3.3 Let I e N, Mq e R-^ , and consider initial data of Section 2.2, 

with M{rj) = {m{r]),rjMQ) and ■= Qj^jf^^y If V is small enough, than the 
outermost apparent horizon is precisely the union of the Schwarzschild horizons 
\x — Xi\ = mi/2. 

Proof: Let Ai be small annular regions around the Xj's, chosen so that the 
metric is exactly Schwarzschild there, then by Theorem 3.2 we have 6"^ C UjAj 
for 77 small enough. The result follows now from the following fact: □ 

Lemma 3.4 The only compact embedded minimal surface in a Riemannian Schwarzschild 
metric (2.3) is the sphere \x — xo\ = m/2. 

Proof: The Riemannian Schwarzschild metric is foliated by spheres of constant 
mean curvature. These are outer mean convex with respect to the normal 
pointing away from the neck. We may now apply the maximum principle. If 
S is any compact embedded (or even immersed) minimal surface, then there 
is some outermost such sphere which makes 'first contact' with S, which is a 
contradiction. The only alternative is that S coincides with one of these spheres, 
and since it is minimal, it must be the central one. 

We may also argue using Lorentzian methods. In fact, standard causality 
theory shows that a compact embedded minimal surface within a time sym- 
metric Cauchy surface cannot be seen from J^+, and so we may obtain the 
conclusion by inspecting the well known conformal diagram for the Kruskal- 
Szekeres extension of the Schwarzschild space-time. □ 
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Using [18, Lemma 4.1] one last time, each component of Si is a sphere, and 
it is plausible that these must agree with the homologically area-minimising 
surfaces Sj C A^, whose topology is a priori either that of a sphere or a torus. 
In each of the examples in the last section, the annular regions Ai are small 
perturbations of rescalings of the Riemannian Schwarzschild metric, and so one 
may construct a foliation by constant mean curvature spheres using the implicit 
function theorem; from this it follows just as before that there is a unique stable 
minimal surface representing cTj, so that S'^ = Ej for all i. However, it is not 
clear that this is true in more general cases. 

There is an analogue of Theorem 3.2 concerning trapped surfaces for asymp- 
totically hyperboloidal initial data sets. Suppose that y has the same topology 
as before, but that the metrics 5^ are asymptotically hyperboloidal. Metrics of 
this sort, with many necks, can be constructed as in Section 2.3. We suppose 
that the diffeomorphism ^'"^ identifies E{r]) with the complement of a finite 
number of balls in (or indeed any asymptotically hyperboloidal manifold 
with constant negative scalar curvature) ; we also replace the hypotheses a) and 
b)by: 

a') [Metric convergence on the distinguished end:] If K is any compact 
subset of \ Ui{xi}i=i,... j, then for some a G (0, 1) 

lim ll*^(5r/) -5^llc2.°(i^) =0; 

here 5^ is the standard hyperbolic metric on H^. 

b') [Neck boundaries with controlled mean curvature:] For in a 

sufficiently small interval (0, the outer boundaries ^i{S'^ x { — 1}) have 
mean curvature h < —2 (with respect to the inward-pointing unit normal). 

We shall be using the maximum principle in the following form. Let and 
S2 be two oriented, connected, embedded surfaces with constant mean curvature 
Hi and H2, respectively. Suppose that these surfaces are tangent at a point p 
and their normals are equal at this point, and that in some small neighborhood 
Si lies on the 'interior' of 5*2 (with respect to the normal). Then necessarily 
Hi > H2, and if Hi = H2, these surfaces must coincide. As a slightly weaker 
statement, if Hi and H2 are now possibly variable and if Hi > H2 everywhere, 
then this one-sided tangency cannot occur. As an immediate application, let 
S be any compact oriented surface in which contains all of the points Xj in 
its interior, and which has mean curvature everywhere greater than —2 with 
respect to its outward normal. (For example, we could let H = S{0, R), a large 
sphere.) This mean curvature remains greater than —2 when computed with 
respect to the metric 5^ when r] is small enough. Hence Srj cannot be internally 
tangent to this sphere, and this shows that in particular 5^ is contained in a 
fixed neighborhood of the convex hull of the x,. 

Proposition 3.5 Under hypotheses a') and b'), there is at least one (smooth, 
embedded, oriented) surface Syj which is homologous to S{0,R) C (for suf- 
ficiently large R) and which has mean curvature —2 with respect to the normal 
pointing into the unbounded component of ^ \ S^, i.e. is marginally trapped. 
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Proof: Since ,y is a manifold with boundary, the volume form dVg^ is exact, 
hence equals dA for some (non-unique) 2-form A. Now define the functional 



Note that changing A alters L by a constant in each homology class, but this 
is irrelevant for our purposes. This functional was studied, for example, in [32], 
and it follows from (2.14) in that paper that if 5 is a smooth stationary point 
of L, then the mean curvature of S is equal to —2. 

Henceforth, let S (0, R) denote any large geodesic sphere in which encloses 
all of the points Xi, and which we identify with a surface in ^ using ^^j. We may 
apply the usual geometric measure theory arguments, as follows, to conclude 
the existence of a smooth minimiser in the homology class of S{0,R). First, 
it is clear that L{S{0,R)) increases without bound as i? — > oo. Next, when 
looking for a minimiser S, we may as well assume that S lies in the bounded 
component U ol ^ \ S{0, R), for if this were not the case, we could replace S 
by a homologous surface S' on which L assumes a smaller value. For example, 
if V is the bounded component oi \S, then d{U PI F) is a suitable^ choice 
for S'. Hence, since wc may assume that any minimising sequence Sj remains 
within a compact set in y, and since L is bounded below, we may find a 
minimiser 5^. The assumption that the outer boundaries have mean curvature 
H < —2 ensures that remains in the interior of S^, cf. [32, Lemma 4]. The 
same regularity theory as was quoted earlier implies that the minimiser Sfj is a 
smooth embedded and oriented surface in the interior of =5^. □ 

Theorem 3.6 Assume g^^) is asymptotically hyperboloidal and satisfies the 
hypotheses a') and h'). For rj in some sufficiently small interval (0, ?7i], any 
trapped surface which is homologous to S{0,R) is contained in Uj^^^Ai and 
has at least I connected components. 

Notice that we are not assuming that iS^ is an outermost trapped surface 
here. 

Proof: We have already indicated that such trapped surfaces exist. To prove 
that Sf] C DiAi, we proceed as before and assume that this is not the case. To 
take a limit as r/ ^ 0, we use the methods and estimates from [24] and [23], 
which adapt in a straightforward way to small metric perturbations of hyper- 
bolic space, cf. also [32, Lemma 2]. In general, the situation is not as simple 
as for stable minimal surfaces because of the possibility of small necks in 5^ 
pinching oif, even in regions where the ambient geometry is uniform. One can 
prove that the limit surface S' is a finite union of smooth embedded surfaces 
Sj which are mutually tangent at their points of intersection. (This part of the 
argument does not use specifically that \H\ =2, and it is possible one could use 

*This follows from convexity: if one lets Si be the portion of S outside the sphere, and 11 
the projection from the exterior onto the surface of the sphere, then n(S'i) has less area than 
Si, because the Jacobian of 11 is everywhere less than 1. So the sphere contribution to L is 
reduced; clearly the volume contribution is reduced as well. 
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this special feature more strongly and show directly that S' is smooth; however, 
this is not so important for our purposes.) We may use the same removable sin- 
gularities theorem as before, or rather its proof, to show that each of the Sj are 
smooth at the points Xj. However, each Sj is compact and has constant mean 
curvature —2. But one could then find a horosphere tangent to Sj, for example 
by bringing it in from infinity (in any direction) until it reaches a point of first 
contact, and this would contradict the maximum principle. Hence Sj could not 
exist. (An alternative nonexistence proof is to note that if such Sj's existed, 
then Minkowski space-time would contain non-empty black hole regions.) 

We have now reduced to the case where S*^ C UAi. The same intersection 
theory argument as in the proof of Theorem 3.2 shows that each of the inter- 
sections Sn n Ai is nonempty, and so Sr) must have at least I components. Note 
that each contains an area-minimising surface Sj which is homologous to 
the outer boundary, and the maximum principle implies that is contained in 
the region between $i(S'^ x {—1} and Sj. □ 

One can impose various geometric conditions on the metric gy on the Ai 
which would ensure that Srj has exactly / components. A rather stringent one, 
which however is satisfied for the asymptotically hyperboloidal initial data sets 
of [19] for 6 small enough, is: 

c') The diffeomorphisms $i can be chosen, now possibly depending on rj, 
so that each sphere $j(S'^ x {t}) has constant mean curvature Hi{t), and 
that each Hi is a monotone function on [—1, 1] with values in some interval 
[—h{r]),h{r])], where h{r]) > 2. 

To see that the initial data sets of §2.3 have CMC foliations on each neck 
region, one can argue as follows. The quantitative estimates for the metric 
grj on these neck regions from [19, §8] show that if we scale {Ai,g^) to have 
a fixed neck size (e.g. to have injectivity radius always equal to 1), then this 
anmilus is quasi-isomctric, with constant tending quickly to 1 as — >■ 0, 
with the neck region for the Ricmannian Schwarzschild space (scaled to have 
the same normalisation). This latter space has a global CMC foliation, and 
by the implicit function theorem we can produce such a CMC foliation in any 
fixed neighborhood of the neck. The outermost leaves of this foliation will have 
mean curvature zth, say, and when rescaled down to the original size, these 
leaves now have mean curvature ±/i(?7), where h{r]) — > oo. 

We use this CMC foliation as follows. Consider the component S'j ,^ — 
Sr, n Ai. Choose r' and r" so that Si^rj C S^ x [t',t"], and such that this is the 
narrowest band with this property. Then Si^r, is tangent to both boundaries, 
and its outward unit normal at these points lies in the same direction as df. 
Denoting by H' and H" the constant mean curvatures of those two boundaries, 
then the maximum principle gives that H' > —2 > H" . But r' < r" and so 
H' = H" and finally Si^^) must coincide with a leaf of the foliation, and hence 
is connected. 
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4 Sections of event horizons have at least / compo- 
nents 



In this section we analyze the global structure of the maximal globally hyper- 
bolic developments of families of initial data sharing certain overall properties 
with those of Section 2.2, when the mass parameters are sufficiently small. This 
question is rather different from the one raised in the previous section, because 
the existence of apparent horizons involves only the geometry of the initial data, 
which is fairly well controlled. On the other hand, the notion of the event hori- 
zon involves the global structure of the resulting space-time, about which only 
very scant information is available. Before proceeding further, the following 
should be said: because gravity is attractive, and because the Schwarzschild 
regions of the initial data of Section 2.2 are initially at rest with respect to each 
other, one expects that those regions will "start moving towards each other", 
leading either to the formation of naked singularities, or to a single black hole. 
In particular the resulting event horizon, if occurring, is expected to be a con- 
nected hypersurface in space-time. Nevertheless, the properties of the maximal 
globally hyperbolic developments {^,g) of the data which we present below 
lead us to conjecture that there exists no slicing of ^ by Cauchy surfaces 5^r 
which are asymptotically fiat in all their asymptotic regions and in which all 
the intersections S"^ fl o,f^ connected. This seems to be the proper way of 
making precise the many-black hole character of certain families of black hole 
space-times. While we do not prove such a conjecture, it follows from what is 
said below that for some configurations there exist natural slicings of ^ which 
do have this property. 

Recall that the black hole event horizon is usually defined as 

S+ ■.= j-{J+;{J,g)) . (4.1) 

Here the causal past J~ is taken with respect to the conformally rescaled space- 
time metric ^ on the completed space-time with boundary ^ := ^ U J^+. 
Thus, the starting point of any black hole considerations is the existence of a 
conformal completion at future null infinity J^'^ . In this context one usually 
assumes that J^"*" satisfies various completeness conditions [16, 17, 31] (compare 
the discussion in [8, 11]). As already mentioned, for the metrics of Section 2.2 
past-completeness of ^+ is guaranteed by the fact that the initial data are ex- 
actly Schwarzschild outside of a compact set. However, the current understand- 
ing of the global properties of solutions of the Cauchy problem for the Einstein 
equations is insufficient to guarantee any future completeness properties of the 
resulting Nevertheless, we shall see that for some of those metrics the 

conformal boundary can be chosen sufficiently large to the future so that 
defined by (1.1) will have more than one component. (This feature will 
persist upon enlarging , and will therefore also hold for a maximal one.) 
Before passing to a proof of this fact let us point out that the existence time 
of the solution, defined as the lowest upper bound on the existence time of all 
geodesies normal to goes to zero as the mass parameters go to zero. In or- 
der to see that, let T be a maximally extended future directed timelike geodesic 
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Figure 1: The maximal globally hyperbolic development of the initial data of 
Section 2.2 in a neighborhood of B{xi,ri). The metric in the dotted region is ex- 
actly the Schwarzschild metric. The 45-degrees sloped dashed lines correspond 
to the Schwarzschild event horizons. The straight part of the boldface lines is 
exactly the Schwarzschild singularity; one expects that some form of that sin- 
gularity will survive in nearby regions influenced by the non-Schwarzschildian 
initial data, as depicted by the curved part of the boldface line, but no results 
of this kind are known. 

normal to ,y starting at the minimal neck of the Einstein-Rosen bridge of the 
usual Kruskal-Szekeres extension {^Schw i dSchw) of the Schwarzschild space- 
time with mass m. Either an explicit calculation, or a simple scaling argument, 
show that the Lorentzian length of T is proportional to m. Now, if S in (2.5) 
is small enough, then the maximal globally hyperbolic development {^,g) of 
the initial data of Section 2.2 will contain a region isometrically diffeomorphic 
to a neighborhood of T in {^Schw, QSchw) as iii Figure 4. This shows that for 
small data g) is necessarily "small" , in the sense made precise above, and 
a complete understanding of the global structure of the resulting space-times 
might be a delicate issue. 

Let us return to the problem of main interest here, namely non-connectedness 
of sections of S^, as defined by (1.1). We shall show that the stability results 
of Priedrich [15] can be used to reduce this question to elementary considera- 
tions of light-cones in Minkowski space-time. Recall that the simplest conformal 
completion of the timelike future of a point in Minkowski space-time (R"^"*"^,//) 
is obtained by performing the space-time inversion 

> , r^apx'^x^ < 0} 9 x'^ ^ y'^ = € {y° < , r?a/3y"/ < 0} . 

(4.2) 

Here rj^u is the Minkowski metric. A drawback of the transformation (4.2) is 
that it does not give the whole conformal completion of Minkowski space-time 
at once; however, a major advantage thereof is that the rescaled metric is again 
the Minkowski one, so that the causal properties of the rescaled space-time are 
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straightforward to analyse, and to visualise: 

Under (4.2) the future timelike cone I~^{Ox; (M^"*"^, r/)) of the origin Ox of the 
coordinates becomes the past timelike cone /"(O^; (M'^+^, ry)) of the origin Oy 
of the coordinates; further, Oj^ is the future timelike infinity point i"*", while 
J~{Oy; {M.^~^^,r))) becomes that part of the Minkowskian ^+ which lies to the 
causal future of 0^; in the conformally completed Minkowski space-time. 

Choose, now, a set of points yi and strictly positive numbers 6i, i = 1, . . . I, 
with 

\yi\ + Si< 1/2 , Si < \yi\ , (4.3) 

with the balls B{yi,di) — pairwise disjoint. The points yi should be thought 
as the y-coordinatcs equivalents of the points of Section 2.2. Let the initial 
surface =5^o be defined by the equation^ 

-^0 = {y° = , < |y| < i , y^B{yi,5i)} (4.4) 

with {yf") = {y^,y). Set 

Ki:={y^ = -^, yGB{yi,Si)}, 

Jir := {-i < 2/° < r} n (M^+i, ^)) \ ( J+(Ki; (M^+i, r?))) , 

:= {-\ < y° < r} n J-(0,; (M^+i,^)) \ J+{Ki- {m?+\r^))) . 

(4.5) 

(See Figures 2-4.) The parameter r should be thought of as the ^''-coordinate 
height of the regions on which the solution associated to the non-trivial initial 
data exists. One can think of the regions J~^{Ki; (M^+^, 77)) as the regions where 
the non-trivial geometry, associated to neighborhoods of the black hole regions, 
is localised. 

In order to be able to take advantage of Friedrich's stability results [15], 
we make the following hypotheses: we consider families of hyperboloidal initial 
data sets {{gT],Kr,)}r,e[0,r,o], defined on B{0, ^)\Ui{yi}i=i^...j, for some ijo e M+, 
/ G N, such that 

a) [Uniform convergence] For any open set which has compact closure 
in B{0, i) \ Ui{yi}i=i^... J we have'^ 

Here gj^ stands for the unit hyperbolic metric. 

^Tho value —1/2 for y'^ is chosen for definiteness; any other value can of course be chosen. It 
is, nevertheless, worthwhile mentioning that this choice corresponds to an upper hyperboloid 
{x° = -1 -I- VI + r^}- It appears that initial data similar to those of Section 2.2 can be 
constructed directly on such hyperboloids by extending the techniques of Corvino and Schoen 
to a hyperboloidal setting. 

®The C'' norms here can be replaced by any Sobolev norms which guarantee that the 
resulting space-time metric, obtained by evolving the initial data using the vacuum Einstein 
equations, is C^. 
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Figure 2: A (2 + 1) -dimensional version of the space-time = 2, for r 

smaller than the time t_ of (4.7). 



b) [Existence of 1 trapped surfaces] There exist r^, i = such 
that for every rj there exists a compact smooth embedded trapped or 
marginally trapped surface Si^r] / satisfying 

Si,n C B{yi,ri) . 

The balls B(yi,ri) will further be required to be pairwise disjoint. In some 
of the arguments below b) will need to be strengthened to: 

b') Moreover, 

limsup{|y - yj| : y £ Sj,.^} = . 

(For marginally trapped surfaces this does follow from a), b), and from 
what is said at the end of Section 3.) 

c) [Existence of The resulting family of space-times admit conformal 
completions which are sufficiently differentiable so that Friedrich's stabil- 
ity theorem (or perhaps some extension thereof, in the spirit of [12,27]) 
applies. 
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Figure 3: The space-time of Figure 2 for < r < r+, with given by 
(4.8); compare Figure 5. 



It is not immediately clear whether the initial data of Section 2.2 are compat- 
ible with those hypotheses. There are a few issues here. Suppose, for instance, 
that the solution associated to the initial data of Section 2.2 remains as close 
as desired to the Minkowski one, when making the mass parameters small, in a 
small neighborhood of the spheres S{xi,ri), for a long time, and that the time 
in question tends to infinity as the mass parameters rrii go to zero. In such 
a case the hypotheses above would obviously hold, whatever the choice of the 
hyperboloidal initial surface =5^o- However, we are not aware of any argument 
which would justify that this is the correct picture, and the discussion around 
Figure 4 suggests that this might actually be wrong. A simple way of avoiding 
the question of the time of existence of the solution near the spheres S'(xj, rj)'s 
is to suppose that all the points Xi lie on the surface of some sphere.^ We 
can then choose the hyperboloid so that the Xj's lie on the intersection of this 
hyperboloid with the hypersurface = 0. For such configurations clearly all 
the hypotheses above are satisfied. 

In this context the following comment is also appropriate: So far we have 

^This involves no loss of generality if / = 2, or if / = 3 and the Xi's are not aligned. 
However, for 1 — 3 the configurations of Section 2.2 are actually co-linear. 
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Figure 4: The space-time of Figure 2 for r larger than the time r+ of (4.9). 



assumed that the initial data are prescribed on the hypersurface given by 
(4.4). The exact choice of .5^q is clearly irrelevant, and a similar picture would 
be obtained at this stage with any hypersurface S^q which asymptotically ap- 
proaches a hypersurface of constant conformal time ip . In particular we could 
choose to coincide with the hypersurface = const for some large posi- 
tive constant for \y\ < R, and to coincide with the hypersurface {y^ = —1/2} for 
\y\ large enough. On =5^o we can use the initial data of Section 2.2 for \y\ < R, 
and appropriate data obtained by time-evolution elsewhere. In such a case, 
the nature of the initial data of Section 2.2 would guarantee that the induced 
data on the new hypersurface would have almost all the properties used in the 
discussion above: The only property missing is that the limiting data set, as rj 
tends to zero, would not be the hyperbolic data set {g = g^, K = g^), but one 
corresponding to an appropriate hypersurface in Minkowski space-time. How- 
ever, our proof of existence of / components of the initial section of the event 
horizons relies on the fact that the radii 6i can be made arbitrarily small on 
a hypersurface of constant y'^-time, and we wouldn't be able to achieve the 
desired conclusions on general hypersurfaces. 

Finally, consider the initial data of Section 2.3 on asymptotically flat hy- 
persurfaces and on hyperboloids; those initial data can be chosen to satisfy 
conditions a), b) and b') above. The initial data of [19] on asymptotically flat 
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hypersurfaces are not known to satisfy condition c). However, those constructed 
in [19] on hypcrboloids can be chosen to satisfy that condition: the resulting 
globahy hyperbolic developments will not have a J^"*" which is complete to the 
past, but it should be clear from the arguments below that this is irrelevant for 
most of the problems discussed here. 

Returning to our model spaces the corresponding model data on ,5^q 
are exactly those for the Minkowski metric. In the physical metric the initial 
data on will be close to the Minkowskian ones, the difference being as 
small as desired when 77 is made sufficiently small. Under conditions a) and 
c) as spelled-out at the beginning of this section, the usual arguments about 
continuous dependence of solutions of hyperbolic PDE's upon initial data over 
compact sets, as applied to the conformal Einstein equations of Priedrich [15], 
show that for any fixed r the physical metric g will exist on Ji^- and will be 
as close as desired to the Minkowski one on ^t-, when the initial data on 

are sufficiently close to the Minkowski ones. This implies that the causal 
structure of the physical space-time on j^-,- will be approximated as accurately 
as desired by that of the Minkowski space-time on when the initial data 
are sufficiently close to the Minkowskian ones on .^q- In particular the figures 
presented here will accurately describe the geometry of null geodesies in the 
physical space-time. 

Let 

^+ := {-\ < y° < r} n j-(0,; (M=^+\r?)) (4.6) 
be the conformal boundary of and suppose that 

r < T_ := -Fminj^ - [yil - 5j| . (4.7) 

In that case the black hole event horizon , in the space-time (^,-,7?), asso- 
ciated with the conformal boundary J"^ ^ 

<r+:= j-(^+;(^„77)) , 

will be a union of spheres: 

= Ute[_ = |y| = -2r + t\ , 

see Figure 5. In particular will be connected, so that each section thereof 
through a hypersurface 

= const} 

will also be connected. This holds for the Minkowski metric, and hence also 
for the physical metric for initial data sufficiently close to Minkowskian ones. 
Thus, if the physical space-time develops a singularity and stops to exist at 
some time r satisfying (4.7), then the boundary of the black hole region will 
be connected. As long as this last possibility occurs it is meaningless — within 
the J'^ framework — to assert that is a multi-black-hole space-time.^ 

*0n the other hand, at an intuitive level it is clear that, whatever the value of r, the 

physical space-time docs contain distinct regions which display "black hole" properties, even 
though this does not fit well into the framework. It seems that any significant insight into 
such situations will be gained only after better understanding of the long time behavior of 
solutions of Einstein equations will have been reached. 
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Figure 5: .-M ^ for t < t^, compare Figure 2. The shaded area is the part of ,y 
which can be seen from J^^^, and its complement in ^ is therefore the (partly 
"physically wrong") black hole region, within S^, with respect to J^^. 



We stress that the global structure of Figure 5 could very well arise for non- 
trivial initial data, whether small or large, even if all singularities are shielded 
by the event horizon (in which case r, near can be thought of as being 
infinite, and should not be identified with a Minkowskian coordinate). The 
point of our considerations below is to show that this will not happen for some 
configurations. 

Now, as soon as the initial value of y*^ exceeds the value r_ given by (4.7), 
some null geodesies starting with this initial value backwards in time from the 
Minkowskian Scri 

^(M3+i,,) := j-(0,;(M2+\r?)) 

enter the region J'^{Ki) where the metric fails to be close to the Minkowski 
one, even for small mass positive parameters mj, and where singularities do 
form in short time. The visibility of those singularities from J^"*" would be 
forbidden if a suitable version of cosmic censorship hypothesis applied, but no 
such results have been established so far. As of today there is no justification 
for the possibility that the physical can be continued uniformly beyond 
the points at which some of the generators of the Minkowskian meet some 
null geodesies emanating from the KiS (though these generators actually do 
continue "a little" in the situation at hand). Whatever the case, stability implies 
that one might continue each generator of the physical J^"*", associated to the 
non-trivial initial data, to the future from the boundary of the initial data 
hypersurface up to the first point the Minkowskian past of which intersects one 
of the KiS. From this point of view the only significant feature distinguishing 
various values of r is that sections of the model J^^, as defined by (4.6), with 
hypersurfaces = t} will be spheres for t < T-, cf. Figures 2 and 5, while this 
will not be the case anymore if r > f > r_. Furthermore, the causal geometry 
of {^rif]) becomes interesting only for 

r > r* := +min|^ - . (4.8) 

We shall not attempt to analyse exhaustively what happens for all t > t^, and 
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Figure 6: with r = r_. The shaded area is the part of 5^ which can be 
seen from J^^. 

all possible values of y^, but we will concentrate on a few specific cases. We set 

#+(0) :=<+n^o . 

We wish to exhibit configurations for which has at least / components. 

Now, by standard causality theory^, the -Sj^^'s of condition b) cannot be seen 
from ^T^, whatever the value of r. It follows that (Si^(O) is never empty. Our 
aim is to construct hypersurfaces jVi C 5^q, i = 1, ...,/ — 1 with the following 
properties: 

\. jVi(Z so that (g;+(o) n ^ = 0. 

2. The jVi'S. separate =5^o into / distinct, open, connected sets Gi such that 
each Gi contains precisely one Sj,^. 

It then clearly follows that <^^(0) has at least I components. 

Let us start with the case 1=2. Without loss of generality one can then 
assume xi = —X2- Further, under the current hypotheses one can without loss 
of generality assume that the constants 5i of (4.3) satisfy 5\ = 82 by replacing 
the smaller of the 5j's by the larger one, and making the parameter q smaller 
if necessary. From now on we assume that rj has been chosen small enough so 
that the physical metric exists on with r larger than 

r+ := - 1 < : (4.9) 

this value of r corresponds to the value of ip at the meeting points of a generator 
of J"'"(Ki, (M'^+-'^, r/)) and a generator of J+(i^2) (IR'^^^, ??)) a generator of 
J~(Oy, (R^"^^,7y)) — see the proof of Proposition 4.1 below. This is also the 
"highest point" of for r > r+, compare Figures 4, 7 and 8. Finally, this 
corresponds to the value of r above which does not change any more: 

Vr > r+ = •^T+ 5 

It is then obvious from Figure 7, in which the K^s are very close to the conformal 
boundary, that the past of the "highest points" of contains points lying on 
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Figure 7: Mq for yi's close to the conformal boundary. 



the straight Une segment connecting the two black holes. This same property 
is still visible, though with a little more effort, from Figure 4 where the black 
hole regions are fairly far away from each other. On the other hand, it should 
be clear from Figure 8 that ^^(0) will be connected there. Before presenting 
a precise form of those statements, let us introduce the notation 



In the proposition that follows the constant —1/2 appearing in (4.4)-(4.5), rep- 
resenting the coordinate of the initial data hypersurface, has been replaced 
by an arbitrary constant tq < 0: 

Proposition 4.1 Suppose that I = 2, 5 := 5i = 82, Vi = -m, set a := 
Then the plane 



^In the case of the initial data of Section 2.2 the use of causal theory is actually not needed: 
the existence of black hole regions follows immediately from the Schwarzschildian character 
of the data on B{xi,ri), as made clear by Figure 5. 



(y^) = {t,x,y,z) . 



^ := {t = To , X = 0} C {t = To} 
is included in J^i-J^^., {^r+jf])) ^^'^ on/y if 




(4.10) 
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Figure 8: Mq for yi's close to each other. 



Remark 4.2 Equation (4-10) always holds for S/\to\ small enough. This is all 
that is needed for our purposes: the sets Ki have been chosen to contain the 
non-trivial geometry, and condition a) guarantees that they can he chosen as 
small as desired by choosing the parameter rj small enough. On the other hand, 
it is restrictive: for example, /or a/|ro| = 1/2, which is the case in Figures 3-4, 
(4.10) leads to 5 < (^/2 - l)|ro|/2, which can fail to he satisfied without violating 
our remaining restrictions 5 < a, 6 + a < \tq\. 

Proof: Clearly pQ := (to,0) G J^{^t+) if and only if the whole plane {t = 
To , = 0} is included in J~ ( J^^^ ) (throughout this proof all the causal objects 
are taken in (M^+^, r/)). Set pi = (tq — (5, y^), then 

J+{Ki) = j+{pi)r\{t>To} . 

Without loss of generality we may assume xi = (a, 0,0). A simple calculation 
gives 

j+(Ki) n j+{K2) n = {(T5,0,y,z) \y^ + z^ = t|} C {x = 0} , 

with 

2(|ro| + 5) 
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On the other hand 

j+{po) n j-{Oy) n{x = 0} = |(^,0,y,z) \y' + z'= (^)'| . 

It then easily follows, e.g. by symmetry arguments, that po will be in the causal 
past of if and only if 

This last equation is equivalent to (4.10). □ 

Proposition 4.1 together with Remark 4.2 settle the case 1 = 2. In order to 
proceed further, it is necessary to understand the geometry of the intersections 

j+{p)nj-{Oy) , pei-{Oy). 

It is convenient to consider general space-time dimensions n + 1. Let tq < 
and let p = (tq, q) G /"(O) C with (f G B{0, \to\) C M"; for the discussion 

here all the causal objects are defined with respect to the Minkowski metric r] 
in R"+-'^. We denote by tt the projection along the first, timelike coordinate axis 
in M"+^ (associated to a coordinate which we denote by x°). We set 

% ■■= J+((ro,gl)n j-(0) , 
:= 7r{%). 

A simple computation shows that the ^g>'s are solid ellipsoids: for q = (a, 0), 
\a\ < I To I, where denotes the origin in M*^"^, we have 

= ^ ' ■= i-r+-H-n' . (4.11) 

Here we use the symbol x to denote the first coordinate in M". For further 
purposes the following properties of the 's are useful: 

• The ^(ao)'^ cigar shaped, except for the one with a = which is 
a ball of radius |to|/2. 

• The ^, ^^'s are centred at (a/2, 0), and their extent in the first coordinate 
X equals tq, independently of a. 

• The intersection of the ff^^ g-j 's with the central hyperplane {x = a/2} is an 

(n — l)-dimensional ball of radius /(a/|ro|), where /(/?) = |ro|i/l — /3^/2. 
The function / : [0, 1] — [0, |to|/2] is strictly decreasing. 

• We further have 

i9^(a,o)) n = a/2} C S'^-^O, |ro|/2) C M" ; (4.12) 

we have decorated the sphere S"'~^{0, |ro|/2) with a subscript n — 1 to 
emphasise its dimension. The sections (4.12) are the "fattest" x-sections 
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of the ^(„o)'^- Thus, as a increases from to |ro| the fattest part of 

^(aO) thinners, with its boundary travehng on 5'"'^^(0, |to|/2) from the 
equatorial hyperplanc {x = 0} all the way to the north pole x = \to\/2. 
The ellipsoids degenerate to a line in this last limit. 

The following simple rules complement the above: 

1. Let X G B^{0, |ro|) \ ^(ao)' Then the point p = {to,x) is in the timelike 
past of 

j-(0)\/+((ro,a,0)) ; 
the required timelike curve is simply a segment of the vertical line t — > 

2. If X G (^'^(a o))'^^(=^t')' then the whole line segment s(a, 0) + (l — s € 
(0, 1) is included in /~(^^): Indeed, by definition of there exists a 

causal geodesic F from (a, 0) to a point p on with p projecting down to 
x — the projection 7r(r) of F is the line segment s{a, 0) + (l— s)x, s G [0, 1]. 
The required timclikc curve is obtained by a timelike deformation of the 
following path: one first moves from (ro,x) towards the future along the 
a;° coordinate line until one meets F, and then one moves along F until 
one meets J^^. We note that if p is not on an edge of J^^, then the whole 
closed segment s(a,0) + (1 — s)x, s G [0, 1], will actually be included in 
/-(^+). 

We consider now the case of three or more Xj's. The simplest configuration 
is that with all the points yi aligned; without loss of generality we can then 
assume that they lie on the axis {y = z = 0}. We have the following result, 
where we do not assume 7 = 3: 

Theorem 4.3 Suppose that all the yi 's are co-linear, and that the hypotheses 
a), h), b') and c) of the beginning of this section hold. Then for rj's small 
enough S'^ (0) has at least I components. 

Proof: The case / = 2 is covered by Proposition 4.1, it is thus sufficient to 
consider J > 3. We will construct hypersurfaces ^ which will be included in 
I~{J^^) for 5i = 0, as all the objects involved depend continuously upon Si the 
result for small (5i's will follow. 

Consider first the case of space-dimension n = 2, let iji = (ai,0). Then 

Suppose that / = 3, without loss of generality we may assume ai > 0, consider 
any point b G (0, ai). What has been said concerning the geometry of the 
^(^aO)^ implies that the set 

% ■■= {^(6,o)n7r(^+)}n{x>o} 

= {^m \ {^(0' l^o|/2) u ^(a„o)}} n {X > 0} 



26 



Figure 9: A (2 + 1) -dimensional version of the space-time = 3, for r 

sufficiently large. 



is non-empty, see Figure 10. Since the dimension is two, the desired hypersur- 
face jV\ is actually a curve, obtained as follows: let x be any point in let 

71 be the segment sx + (1 — s)(6, 0), s G [0, 1]. Let 72 be the path obtained by 
first following 71, and then a line parallel to the y axis. By the rules 1) and 2) 

72 is included in /^(.^^+ ). We define JV\ to be the union of 72 and of its image 
under the map — > {x,—y). The second hypersurface ^ is obtained by 
taking the image of ^yi^i under the map {x,y) — > {—x,y). 

It should be obvious to the reader how this construction generalises to higher 
/'s. We simply note that for / = 2N one of the curves, say can always be 
chosen to be the axis {x = 0}. Figure 11 illustrates the case 1 = 4. 

For dimensions n > 3 the desired hypersurfaces can be obtained by rotating 
the curves constructed above around the axis {y^ = • • • = = 0} using the 
action of 50 (n - 1) C SO{n). □ 

Using our discussion of the geometry of the i^q's, with a little work the 
reader should be able to verify the following: 

Theorem 4.4 Let |nj| = 1 for all i, and let yi = An^. There exists < Aq < 1/2 
such that for all Aq < A < 1/2 the black hole region (^^(0) has at least I 
components. 
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Figure 10: Three black holes, aligned. The outermost circle represents the 
conformal boundary of o?o- The dotted region is the projection on the initial 
hypersurface of the part of J{^y) which has been excised by the removal of 
J+(Ki) U J+(if2)- The shaded region is fb. 




Figure 11: Four black holes, aligned. The dotted region is the projection on 
the initial hypersurface of the part of J{Qy) which has been excised by the 
removal of J+(Ki) U J+(i^2) U J^K^). 

□ 

One expects that there exist configurations with / > 3 for which ^^^(O) 
has less than / components. While it is easy to imagine such configurations, a 
justification does not seem to be straightforward. 

So far we have not been assuming anything about how long the metric 
remains close to the fiat one in a small neighborhood of the spheres S{xi,ri). 
The longer this happens, the larger the set of slices {?/° = r} at which the 
horizon has more than one components. This is made clear by Figures 12-13. 
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